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DETERMINATION OF WORK FUNCTIONS NEAR MELTING POINTS OF 
REFRACTORY METALS BY USING A DIRECT -CURRENT ARC 
by William A. Gordon and Gilbert B. Chapman II 
Lewis Research Center 

SUMMARY 

A direct-current arc in argon at atmospheric pressure was used to determine effec- 
tive work functions of refractory metals, including tantalum , tungsten, molybdenum, and 
niobium. The procedure is experimentally advantageous because surface cleanliness of 
the specimen is not critical, high vacuum is not required, and the anode-cathode spac- 
ing is not critical. The experimental procedure involves striking an arc to a metal wire 
cathode to form a melted ball having an emitting area defined by its diameter. The liter- 
ature melting point of the metal is taken as the emitting temperature. By using these 
parameters and the known arc current, effective work functions were calculated from the 
Richardson -Dushman equation. The calculated work functions agree with recommended 
handbook values to within about 0. 1 volt and have typical repeatabilities of 0. 02 volt. By 
varying the arc current, Richardson plots can be made over a temperature range from a 
few hundred degrees below the melting point to about 50° over the melting point of the 
test metal. A Richardson plot over this temperature range is presented for tantalum. 


INTRODUCTION 

The free-running direct -current arc at near atmospheric pressure is generally con- 
sidered to be characterized by thermionic emission of electrons from the cathode (ref. 1). 
However, the arc has not been used as a method for quantitatively determining therm- 
ionic work functions. One literature reference was found which alluded to this possibility 
(ref. 1), but quantitative calculations were not reported. Current reviews of various 
methods of determining work functions do not refer to the direct -current arc (refs. 2 
and 3). 

Recently, the first author developed an argon arc configuration for use in spectro- 
chemical analysis (ref. 4). Conditions were reported for achieving spatial stability of 



the arc by means of a cathode with a spherical tantalum tip. In the development of this 
cathode, studies were made of various refractory -metal cathodes with respect to cathode 
temperature, arc current, and cathode emitting area. Under certain conditions, the re- 
lations among these parameters were reproducible and consistent with the Richardson - 
Dushman equation 


1 " V 2 exp (^) (1) 

where 

2 

j emitting current density, A/cm 

2 2 

A q Richardson emission constant (120.4 A/(cm )(K )) 

T cathode temperature, K 

IQ 

e elementary charge (1. 602x10” C) 
cp work function, V 

k Boltzmann constant (1. 380x10”^ J/K) 

Further investigation has demonstrated that the arc discharge can be advantageous 
for determining thermionic work functions of refractory metals at, or near, their melt- 
ing points. This method is simple and a measurement can be made in about 5 minutes. 
Some other advantages of the arc method are negligible contamination of the emitting sur- 
faces, high vacuum is not required, and electrode spacing of about 1 centimeter is not 
critical. The determinations are highly repeatable, and the work functions determined 
agree with the recommended literature values. 

The results of this work are explained by using a simple model of the arc which as- 
sumes that the electron emission from the cathode is purely thermionic and that all the 
arc current is electronic. Possible errors due to positive ion currents and electric 
fields are also discussed. 


APPARATUS AND MATERIALS 

The apparatus used in this work consisted of an airtight chamber, a direct -current 
arc supply, refractory-metal wires, and an optical pyrometer. The relevant character- 
istics of these items are described in this section. 
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Airtight Chamber 


The airtight chamber is shown schematically in figure 1. It consisted of a quartz 
cylinder sandwiched between steel plates with provisions for introducing gases and for 
evacuation. An airtight seal was made between the quartz ring and the steel plates by 
means of rubber gaskets recessed in grooves machined in the steel plates. The plates 
were connected to the electrical supply with the test wire as the cathode, as shown. A 
bellows seal in the cathode retractor allowed the spacing between the electrodes to be var- 
ied from contact to about 15 millimeters. 



Figure 1. - Schematic of arc chamber for work function determinations. 


Direct -Current Arc Supply 

The arc supply was full -wave rectified, thyratron controlled, and inductor -capacitor 
filtered. The root -mean-square ripple was about 3 percent at 30 amperes. The primary 
voltage was 220 volts at 60 hertz. The secondary voltage applied to the electrodes was 
about 180 volts. Currents from a few amperes to about 60 amperes could be supplied. 
The current to the arc supply was controlled by a feedback loop to ensure constant cur- 
rent output irrespective of electrode spacing or vapor composition in the arc column. 
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Refractory-Metal Wires 


The polycrystalline metal wires used in this work and their diameters in millimeters 
were as follows: tungsten (W), 0.3; rhenium (Re), 1.75; tantalum (Ta), 0.75; molyb- 
denum (Mo), 1. 0; and niobium (Nb), 1. 0. These materials were of ordinary metallur- 
gical grade and were nominally between 99. 5 and 99. 9 percent pure. 


Optical Pyrometer 

Temperature measurements for the data used to make Richardson plots were made 
with a micro-optical pyrometer located about 1. 52 meters (5 ft) from the arc. The 
method of calibrating the pyrometer and correcting for the optical effects of the quartz 
cylinder and the emissivity of the cathode metal are cited in the section PROCEDURE. 


PROCEDURE 

Three variations of the procedure for determining work functions are described in 
this section. These procedures were used to determine work functions from (1) current 
density measurements at the cathode melting point; (2) a plot of the linear form of equa- 
tion (1), commonly called a Richardson plot; and (3) the slopes of Richardson line seg- 
ments as a function of temperature by a more detailed analysis of the data from proce- 
dure 2. Method 1 was applied to Nb, Mo, Ta, W, and Re wires. Methods 2 and 3 were 
applied to Ta wires only to illustrate the procedure by which the effective thermionic 
emission constant, as well as the work function, can be determined. 


Determination of Work Functions at Cathode Melting Point 

The refractory -metal wire to be tested was mounted as the cathode, in contact with 
an anode configuration, as shown in figure 1. The anode was a piece of sheet stock, of 
the same metal as the cathode, resting on a brass block. This anode configuration pro- 
vided conductive cooling, reducing the evaporation of the anode material. 

3 

The chamber was evacuated to about 10 torr and backfilled with argon to about 700 
torr. The arc current source was adjusted to deliver about 30 amperes, the precise cur- 
rent being read from the current meter during operations. The arc voltage was applied; 
and, after a few seconds with the electrodes in contact, the wire cathode was withdrawn 
to form the arc gap. The final arc gap was adjusted to about 1 centimeter. The cathode 
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spot formed on the wire tip and fused the wire to form an approximately spherical emit- 
ting area. The arcing was continued until the sphere ceased to increase in size. Nor- 
mally, about 30 seconds of arcing was sufficient to achieve the maximum sphere size. 
When no increase in sphere size was observed over a period of 15 to 20 seconds, the 
sphere size was assumed to be at thermal steady state. Once this had occurred, the 
emitting sphere would not further increase in size. The arc was then terminated, the 
wire cathode was removed from the arc chamber, and the diameter of the spherical tip 
was measured. 

The work function calculations were made by using the Richardson-Dushman equa- 
tion (1). The current density j was determined by the ratio i/a, where i is the arc 

current and a the emitting area as determined from the measured sphere diameter by 
o 

the relation a = 7rd . The emitting temperature T is the metling point of the cathode 
metal, taken from reference 5. 


Richardson Plots 


Equation ( 1) can be written 


ln Q=-<5) + InA (2) 

where A is the effective emission constant as distinct from A„, the fundamental emis- 

2 ° 

sion constant. Therefore, by plotting In (j/T ) as a function of e/kT, a linear graph 
was obtained with a slope of negative cp and an intercept of In A. 

Richardson plots require data of emitting temperature T as a function of emitting 
current density j. Such data were obtained by first forming a molten spherical emitting 
surface at the tip of the refractory -metal test wire while running the arc at a moderately 
high current, for example, 50 amperes for tantalum. Then, the arc current was reduced 
at intervals of about 2 amperes while the cathode emitting temperature was read at each 
current interval with a micro-optical pyrometer. After reaching a lower current of about 
10 amperes, the arc current was then increased and the emitting temperature again re- 
corded at about 2 -ampere intervals up to about the starting current. 

The temperature readings were obtained with a micro -optical pyrometer which was 
calibrated by using a National Bureau of Standards tungsten ribbon lamp placed inside the 
quartz cylinder to correct for optical effects. The pyrometric temperature readings 
were then corrected for the emissivity of the tantalum (0. 36) from reference 6, using a 
method of correction from reference 7. 

This procedure produced the emitting temperature and current data, at constant cath- 
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ode emitting area, over a temperature range from about 3000 to 3300 K. These data were 
treated first by computing a least-squares linear fit having a slope of -cp and an inter- 
cept equal to In A. A quadratic fit was also determined to detect significant deviation 
from the predicted linearity. Next, a FITLOS (ref. 8) curve -fitting procedure was ap- 
plied to examine in more detail what appeared to be nonlinear distribution of data at the 
high-temperature end of the plot. The force-fitting provision of this program was used, 
allowing a 1 -percent deviation of the data from the fit for each segment of the curve. The 
necessary number of segments and locations of spline joints were determined by the com- 
puter’s utilization of the data. This procedure allowed detection of significant deviations 
from the least -squares linear and single quadratic fits in various segments of the total 
range of data. 

Such a fit, being a smooth coupling of up to 10 quadratics, can follow the data more 
closely through its nonlinear excursions. It thereby is a fairly good detector of regions 
of nonlinearity. Deviations of this fit from linearity were indicated by significant varia- 
tions in its first derivative. 

This experiment was conducted three times, and the data obtained were processed in 
the manner described in this section. 


RESULTS 

Effective Work Functions 

Effective work functions were determined for metals including W, Ta, Mo, Re, and 
Nb. Table I is a summary of effective work functions obtained from measurement of 


TABLE I. - SUMMARY OF WORK FUNCTION 
DETERMINATIONS BY DIRECT -CURRENT 


ARC METHOD 


Metal 

Melting 
point, a 
K 

Work functions, V 

By arc method 

Recommended value* 3 

Mo 

2883 

CO 

4.3 

Nb 

2741 

4.0 

3.99 

Re 

3453 

5. 1 

5.0 

Ta 

3269 

c 4. 13 

4. 12 

W 

3683 

4. 5 

4.5 


a Ref. 5. 
b Ref. 2. 

c Standard deviation for five consecutive determina- 


tions, 0. 013 V. 


6 




sphere diameters formed at 30 amperes and by application of equation (1). Also listed in 
table I are the recommended values from reference 2 for these metals. Metal cathodes 
formed in typical experiments are shown in figure 2. The repeatability of the work func- 
tions for Ta was 0. 013 electron volt, standard deviation, for five consecutive determina- 
tions. This repeatability was typical of the other refractory metals. The work functions 
were also in good agreement with the recommended values for refractory metals from 
reference 2. 



w Ta Re 

Current. 34 A Current, 30 A Current, 10 A 

Diameter, 0.864 mm Diameter, 1.57 mm Diameter, 2. 11 mm 



Figure 2. - Examples of spherical emitting surfaces formed on refractory-metal cathodes. 
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Determination of <p and A by Richardson Plots 


Richardson plots, from which cp and A values for equation (1) can be determined, 
can also be made by the arc method. Data for tantalum are shown in figure 3 with a lin- 
ear least -squares fit (solid line) and a quadratic fit (dashed line) superimposed. The 
standard deviation of the quadratic fit (0. 085) is not significantly better than that of the 
linear fit (0. 093), thus indicating the linear fit to be as precise as the higher degree fit. 
This linear least-squares fit was used to determine and A in the Richardson- 

p 

Dushman equation (1), which predicts a linear relationship between In j/T and e/kT, 

having a slope of -cp and an intercept of In A, in equation (2). From such a fit of the 

data in figure 3, over the full temperature range shown, <p = 4.50 volts, and A = 343. 9 
o o 

A/(cm )(K ). Another similar experiment was conducted; and when the data obtained 
were treated in the manner described, yielded <p = 4.46 volts and A = 311. 04 
A/(cm 2 )(K 2 ). 

Applying a FITLOS curve -fitting procedure to the data resulted in the curve having a 
standard deviation of 0. 082 shown in figure 4(a). This sectioned fit revealed a region of 



3300 3250 3200 3150 3100 3050 3000 

Temperature, T, K 


Figure 3. - Richardson plot for tantalum near melting point. Work function 9 , 4.50 volts 
(slope = - 9 ); effective emission constant A, 343.9 A/(cm 2 )(K 2 ). 
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maximum inflection in an otherwise smooth, near -linear fit. Another presentation, show- 
ing the first derivative of the fit to these same data, appears in figure 4(b). Here the sig- 
nificance of the inflection point, relative to the average derivative, and its precise loca- 
tion on the abscissa are clearly demonstrated. The negative slope cp and the anti -log of 
the intercept A obtained from a linear least-squares fit of the data on the lower temper- 

9 o 

ature side of the discontinuity were 4.04 volts and 60.64 A/(cm )(K ), respectively; and 

9 

a similar fit on the higher temperature side yielded (p = 4. 80 volts and A = 1072 A/(cm^) 
(K^). The fits of both segments have standard deviations of 0. 04. Data obtained from two 
other similar experiments yielded results comparable to the data of figure 4. 



i i i i i i i 

3300 3200 3100 3000 

Temperature, T, K 


(b) First derivatives of fits. 

Figure 4. - FITL0S smoothing of Richardson plot and 
sectioned linear fit. 
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DISCUSSION OF RESULTS 


The work function determinations made at the melting point of the cathode materials 
exhibit good repeatability and apparent accuracy. An important factor in achieving good 
repeatability by this method is the use of a precisely fixed emitting temperature, which 
allows use of the well-known melting point. An examination of equation (1) shows that the 
emitting temperature is the most critical variable because it appears in the exponential 
term. Another factor in the experiment which is conducive to good repeatability is the 
formation of a melted emitting sphere, which provides a well-defined emitting area. 

The experimental procedure for determining work functions by the arc method is rel- 
atively simple. Only nominal evacuation of the airtight chamber is required, and no 
special preparation of the specimen surface is necessary. Other experimental param- 
eters, including the anode -cathode spacing and the argon pressure in the chamber, are 
also not critical. 

The interpretation of these work functions measured from liquid metal surfaces must 
be considered in light of theories for the solid state. It is well known that work functions 
are dependent on the crystalline structure of the metal surface. For melted cathode sur- 
faces obtained in this work, the work functions can be considered to be representative of 
an equivalent solid crystalline form, for example, a hexagonal close -packed structure 
(ref. 9). Following this reasoning, the work functions calculated for liquid surfaces of a 
refractory metal such as tantalum can be considered to be equivalent to the (110) plane of 
a body-centered cubic metal. However, no calculations were made along these lines; and 
further discussion of the theoretical interpretation of these data is outside the scope of 
this report. These remarks, however, are intended to emphasize a problem in interpret- 
ing the results reported herein and suggest a possible approach to explain these results 
for liquid surfaces in terms of the solid state. 

In spite of the problems of theoretical interpretation of these data, the consistency 
and apparent accuracy of the effective work functions suggests that this technique can be 
useful for investigating thermionic parameters of metals and alloys at very high 
temperatures. 


Work Functions and Effective Emission Constants from Richardson Plots 

A comparison of the linear and quadratic graphs shown in figure 3 , and of their fit 
precisions, showed only minor differences, indicating that the data essentially conformed 
to equation (2). The data used to plot the Richardson line are subject to systematic pyro- 
metric errors. Although the pyrometer calibration error was less than 20 K, including 
the correction for the quartz cylinder, the true emissivity of the curved Ta surface at 
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3000 K is now known. The value of emissivity from the literature (0. 36) used to correct 
the pyrometer temperature readings only approximates the emissivity under the exper- 
imental conditions. However, from comparison of the experimental temperature of the 
melted sphere (3230 K) with the literature value for the melting point of Ta (3269 K), we 
estimated a probable total temperature error of about 50 K. 

Examination of the FITLOS plot in figure 4 revealed that the most distinct deviation 
from linearity occurred within a relatively small region of e/kT. The sharp variation of 
the first derivative in that region and the nearly constant derivatives on either side of this 
maximum inflection point are indicative of an abrupt change in the thermionic emission 
from an otherwise well-behaved thermionically emitting surface. The truest fit, in terms 
of the fit precision, was obtained when the data were fit in two linear segments - one 
above this inflection point and one below. 

The interpretation of the inflection point in the Richardson plot is believed to be re- 
lated to a solid -liquid transition of the emitting surface as a function of arc current. 

When the cathode tip was first formed by melting, the interface between the solid wire and 
the melted sphere was approximately at the melting point of the cathode material. How- 
ever, the average sphere temperature was slightly higher because of the small temper- 
ature gradient existing between the sphere surface, wherein the arc energy is input, and 
the wire -sphere interface. Therefore, when the arc current was reduced a few amperes 
below the initial current used to form the sphere, the emitting surface froze, resulting in 
an abrupt change in emissivity and in apparent emitting temperature. The difference be- 
tween the melting point and the average sphere temperature was about 55 K, as indicated 
by the temperature range of the upper curve segment in figure 4(a). 

The small temperature gradient associated with the emitting surface was the basis 
for interpreting the inflection point of the Richardson plots. The change in emissivity ac- 
companying the freezing or melting of the emitting sphere caused a shift in the apparent 
temperature as measured with an optical pyrometer. This systematic change in emissiv- 
ity qualitatively accounts for the displacement of the two sections on the graph. The work 
function determined from the slope of the lower temperature section (higher e/kT) of fig- 
ure 4(a) is not necessarily the same as those determined by calculation from equation (1). 
As pointed out in reference 2, the slope of the Richardson line is a measure of cp Q , the 
work function at absolute zero, and does not reveal a temperature dependency dcp/dT of 
the work function. Such a dependency is usually constant and is, therefore, indistinguish- 
able from a change in the A value. The slope and intercept of the higher temperature 
section of the Richardson plot of figure 4(a) are significantly greater than the slope and 
intercept of the lower temperature section, apparently indicating a change in d cp/dT of 
Ta at the melting point. This explanation is in agreement with reference 10, in which 
both experimental results and a theoretical derivation are given for copper to substantiate 
the existence of a discontinuity of dcp/dT at the melting point. In our experiment, visual 


11 



observations also tended to support this interpretation. The appearance of the emitting 
surface became visually duller as the current density was lowered through the level corre- 
sponding to the observed maximum inflection point on the Richardson plot. Nevertheless, 
it was not possible in the experiment to directly detect the solid -liquid transition instru- 
mentally. Therefore, direct experimental evidence establishing unambiguously that the 
inflection point of the Richardson plot was indeed a result of a phase transition could not 
be obtained. 


Thermionic Model of Direct-Current Arc 

The calculation of work functions by the arc method will be accurate only if the cath- 
ode emission in the direct -current arc is quantitatively described by the parameters of 
equation (1). Insofar as the cathode emission might deviate from this simple model, the 
work functions can be inaccurate. Although proof of the accuracy of the work function 
measurements is not within the scope of this report, it is instructive to qualitatively ex- 
amine possible sources of errors and, where possible, to estimate their magnitude on 
work function determinations. Because of the simplifying assumptions made in the 
method of calculating work functions described herein, it might appear surprising that the 
arc technique has even an apparent validity. The straightforward calculations of work 
functions from the Richardson-Dushman equation do not take into account energy paths 
such as positive ion currents and electric fields. However, on closer examination, it is 
possible to postulate a model of the arc in which the accumulative error in work functions 
is less than about 0. 1 electron volt. In discussing the effects of various error sources, 
it is necessary to describe the arc in more detail. 


Arc Characteristics 

The arc used in this work was a free-running and self -sustained discharge in a static 
argon atmosphere with nominal purity of 99. 99 percent. The pressure was about 700 torr 
but was not critical. The cathode -to -anode distance was about 12 millimeters and the gap 
voltage was about 14 volts. The classification of this arc on the basis of voltage -current 
characteristics compared with other arc types is shown in figure 5. Figure 5(a) is a 
typified diagram showing the relative location of the arc used in this work with respect to 
other gas discharges. The main curve, in figure 5(b), is the static voltage -current rela- 
tion of this arc. Arcs of this type, in which the voltage is relatively independent of cur- 
rent, are typically collision -dominated discharges having ion and electron densities of 
14 15 

the order of 10 to 10 per cubic centimeter (ref. 11). However, the arc column com- 
position is inhomogeneous, especially in the region of the cathode. In this region there 
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(b) Experimental characteristics of argon arc used. 

Figure 5. - Voltage-current characteristics of gaseous discharges. 



Aril Aril Aril 

349. 154 nm 354. 584 nm 358. 844 nm 


Figure 6. - Spectra with arc column focused on spectrographic slit showing argon ion lines excited near cathode. Pressure, 700 torr (120 kN/m^i; 
magnification, approximately 1. 


is an ion sheath, shown in the photographed spectrum of figure 6. The ion space charge 
can promote electron emission under saturation conditions and may also cause errors in 
work function measurements through formation of local electric fields. Furthermore, 
the presence of an ion sheath raises the possibility of errors due to positive ion current 
because in the use of equation (1) it was assumed that all the current is carried by therm- 
ionic electrons. 

Because electrical neutrality of the plasma requires the presence of positive ions, 
the assumption that the current is totally electronic cannot be rigorously valid. In the 
following sections, the inaccuracies of work function determinations caused by deviations 
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from the simple arc model are discussed. From data available on similar arcs it can be 
tentatively shown that the error accruing from use of the simplified model is less than 
about 0. 1 electron volt. 


Electric Fields (Schottky Effect) 

The presence of electric fields in the vicinity of the emitting surface can lower the 

measured work function by an amount equivalent to the work done on the electrons by the 

field. The Schottky effect describes the reduction of the work functions by the field and 

can be expressed by Acp = 3. 79x10 E^'^, where Acp is the change in apparent work 

function caused by the field E in volts per centimeter (ref. 3). 

As can be seen by examination of this expression, to exert a significant influence on 

the total, the electric fields must be higher than could conceivably exist in an arc with a 

4 

total voltage drop of only 20 volts per centimeter. In fact, a field of 10 volts per centi- 
meter would reduce the effective work functions by only 0. 038 electron volt. Therefore, 
the Schottky effect can be dismissed as a significant source of error in work function de- 
terminations by the arc method. 


Positive Ion Currents 

As mentioned previously, the contribution of positive ion currents at the cathode 
surface must be subtracted from the total arc current in order to obtain a proper value 
for thermal electron current for equation (1). Although it is widely accepted that elec- 
trons carry the major portion of the current in direct -current arcs, it is also certain 
that some portion of the current is carried by positive ions. This derives from the fact 
that the primary heat input to the cathode is by ions falling through the cathode field and 
impacting on the cathode surface. To ignore positive ion currents, as was done in this 
work, can result in an erroneously low work function. To estimate the magnitude of pos- 
itive ion currents and the resulting error in work functions, it is necessary to quantify 
the energy balance at the cathode surface. This is a problem of such practical and theo- 
retical difficulty that no method has been devised to obtain this information in arcs with 
thermally emitting cathodes. In this section we shall discuss two models of the cathode 
energy balance from which estimates of positive ion currents can be made. 

Voltage -dominated model. - Estimates of positive ion currents have been made by 
using a voltage -dominated model of the arc. Here we will consider how these estimates 
are made in order to place limits of error on the work functions determined by the arc 
method. A discussion of the approaches used in making these estimates is summarized 
in references 11 and 12. Estimates of positive ion currents in arcs by various workers 
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range from 1 to 36 percent of the total current. A typical treatment of energy relations 
writes the energy input to the cathode as ip(V c + V.), where i^ is the positive ion cur- 
rent, V is the voltage drop in the cathode fall, and V- is the ionization potential of the 
gas atoms. The heat absorbed on recombinations of ions in the impact process can be 
written ip<p, where cp is the electron work function of the cathode surface. Therefore, 
the net energy input is i (V + V. - <p ), whereas this must be balanced by electron cool- 
ing expressed by i & <p, that is, the product of the electron current and the electron work 
function. The energy balance is thus approximated by i <p = i (V + V, - cp). The pos- 
itive ion current fraction can be written ip/(ip + i e ) = <p/(V c + Vj). 

The value of V c can be estimated as about one -half of the total arc voltage, or about 
10 volts. The ionization potential of argon Vj is about 16 volts, and cp for refractory 
metals is about 4 volts. Using these values, the positive ion current can be estimated at 
about 20 percent of the total current. If 20 percent is taken as an estimate of positive ion 
current, the work functions derived by the arc method are low by about 0. 06 electron 
volt. 

As discussed in reference 11, this treatment neglects some potentially significant 
energy paths - notably loss of ion energy by collisions between atoms in the cathode fall 
region, and kinetic energy loss by recoil of neutralized ions at the cathode surface. If 
these processes are important, the work function errors are larger than estimated here- 
in. On the other hand, additional heat input to the cathode can be derived by thermal con- 
duction of the hot plasma to the cathode surface. Because of this effect the positive ion 
current has been estimated to be as low as 2 percent (ref. 13), a negligible amount for 
the purpose of this work. The energy exchange model in which the thermal energy pre- 
dominates is discussed in the next subsection. 

Collision-dominated model. - In the preceding analysis of positive ion currents, it 
was assumed that the voltage fields near the cathode surface provided the driving force 
for the energy distribution among electrons, atoms, and ions. However, at the high tem- 
peratures and high particle densities existing near the cathode, the thermal motions of 
the particles must be considered. It is rather well established for arcs of the type used 
in this work, that random thermal motions of the particles are predominate over voltage 
fields in determining arc mechanisms. This is to say that, as a first approximation, the 
electrons, atoms, and ions are in local thermodynamic equilibrium (LTE) (ref. 14). If 
LTE also exists within a few mean free paths of the cathode surface, the relative propor- 
tion of positive ion current to the total current i^/(i + i g ) can be estimated from the con- 
servation of momentum and energy as m^/^p^ + m*/^. IPor ar & on i° ns this ra ti° 
is very small, indicating negligible positive ion currents and, therefore, negligible error 
in calculating work functions by the Richardson-Dushman equation. 

Although this model can be supported experimentally for macrovolumes in the arc, 
it is not necessarily valid in microvolumes near the cathode. The existing data on ther- 
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mal populations were obtained spectroscopically, and it is not feasible to make such ob- 
servations within a few mean free paths of the cathode surface. Without such direct ex- 
perimental observation, the question of errors due to positive ion currents remains un- 
resolved. However, from the foregoing discussion, when a voltage -dominated model is 
used, 0.06 electron volt would appear to be a reasonable estimate of the maximum error 
in the work functions calculated herein. 


Observations in Helium Atmospheres 

All the data presented in this report were obtained in argon atmospheres. However, 
some experiments were conducted in helium atmospheres at higher pressures to deter- 
mine the criticality of experiment conditions on the experiment results. Substituting he- 

o 

lium for argon and raising the gas pressure to about 4500 torr (600 kN/m ) resulted in 
rather drastic variations of discharge conditions. The appearance of the arc in helium 
is quite different than in argon - the most notable effect being the shifting of the ion 
sheath away from the cathode region to near the anode. In addition, more energy is lost 
by thermal conduction in helium because the heat conductivity is about 10 times that of 
argon. It is also well known that the thermophysical properties of the arc are radically 
different in helium as compared with the arc in argon because of the higher excitation and 
ionization potentials of helium. Despite changes in discharge conditions caused by helium 
atmospheres, and by the higher pressures, the work functions for tantalum in helium 
were not significantly different than those obtained in argon. 


CONCLUDING REMARKS 

The use of a direct -current arc in an argon atmosphere provided a simple means for 
determining work functions of refractory metals in the temperature range near their 
melting points. The work functions are in agreement with recommended handbook values 
to within about 0. 05 volt. The repeatability of the methods resulted from the operation 
of the cathode at precisely its melting point, which allowed use of melting point data for 
the emission temperature in the Richardson -Dushman relation, and from the formation 
of a well-defined emitting area. With this procedure it was also possible to determine 
work functions and effective emission constants by means of Richardson plots. These 
plots exhibit an inflection point at temperatures corresponding to the cathode melting 
points. 

A need for further work on this procedure is indicated. The work functions deter- 
mined on liquid surfaces need to be interpreted in terms of an equivalent crystalline 
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orientation in the solid state. In addition, a more quantitative analysis of energy trans- 
fer mechanisms in the arc region near the cathode is necessary to independently estab- 
lish the accuracy of the method. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 8, 1972, 

134-03. 
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